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' Abstract. We consider portfolio optimization in futures markets. We model the entire futures price curve 

' at once as a solution of a stochastic partial differential equation. The agents objective is to maximize her 

^ ' utility from the final wealth when investing in futures contracts. We study a class of futures price curve 

' models which admit a finite-dimensional realization. Using this, we recast the portfolio optimization problem 

^ as a finite-dimensional control problem and study its solvability. 

^ ■ 1. Introduction 

Cl ; _ . 

' : ^ Futures contracts, see, e.g., [8], form an important class of financial instruments exchanged on commodity 

? markets. These contracts convey the right to purchase or sell a specified quantity of the commodity at a 

O^' fixed price on a fixed future date. This right has to be exercised if the contract is held until the maturity. 



The price fixed in contract is called the futures price and it is set such that no money is paid upfront, i.e., 
the initial value of a futures contract is zero. However, a futures contract yields a cash flow during its life 
time generated by the changes in futures price over time. More precisely, the party in whose favor the futures 
' price change occurs must immediately be paid the full amount of the change by the losing party. Typically, 
r a margin account is set up for this purpose. 

The objective of this paper is to study portfolio optimization on commodity futures markets using tools 
from mathematical finance and stochastic calculus. We assume that the futures prices are settled continuously 
in time and that there is a liquid market for futures contracts for every time-to-delivery y > 0. Now we are 
facing a portfolio optimization problem in infinite dimensions, since we have a stochastic variable, that is, 
a futures price, for a continuum of times to delivery. To tackle this problem, we start with reparametrized 
futures price dynamics a-la Musiela which allows to formulate the dynamics as a solution of a stochastic partial 
differential equation. This equation can be regarded as a stochastic differential equation taking values in a 
real separable Hilbert space. After this, we focus on futures prices models which admit a finite-dimensional 
realization. Roughly speaking, for a given instantaneous futures price model t H> /t, we look for a sufficiently 
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regular function (j> : R+ x R and a c?-dimensional diffusion Z such that ft{y) = 0(2/, Zt) for all t > 0. 

Here, d is the dimension of the realization. Various classes of models admitting finite-dimensional realization 
have been studied over the recent years, mostly in connection with term structure models of interest rates. In 
a series of papers including [3l [H [S] |6] the authors study geometric aspects and finite-dimensional realizations 
of interest rate models using differential geometry and systems and control theoretic methods. In another 
series of papers including [T^l [S], see also [TT], the authors analyze the properties of finite-dimensional 
realizations using the so-called convenient analysis on Frechet space. For example, they prove that the 
function (j) introduced above is necessary affine. In these papers, they provide a general theory on invariance 
of manifolds for infinite-dimensional stochastic equations and study the term structure models of interest 
rates as an application. While being mathematically very general, their analysis is quite demanding to the 
reader. Recently, a more direct approach to affine realizations of term structure models was introduced in 
[22] - we follow this approach when studying our futures price model. We also refer to [10, and [2T|, where 
portfolio optimization with an infinite dimensional state variable is studied. These studies are concerned with 
optimal bond portfolios where the underlying interest rate term structure follows a solution of a stochastic 
partial differential equation. In addition to addressing a whole different application, these studies operate on 
the infinite dimensional level, whereas we study models with a finite-dimensional realization. 

Our paper makes two main contributions. First, we provide a general mathematical framework to treat 
futures portfolio optimization using the up-to-date theory of the realizations of stochastic partial differential 
equations. More precisely, we provide conditions under which a given infinite-dimensional portfolio optimiza- 
tion problem can be solved in terms of a finite-dimensional control problem. Here, the finite-dimensional 
realization (more precisely, the driving coordinate process Z) of the futures price curve t t-^ ft plays a key 
role. Furthermore, we discuss an economic interpretation of the coordinate process and how a solution of 
the finite-dimensional control problem can be connected to the coodinate process and, consequently, back to 
the infinite-dimensional portfolio problem. Second, we extend some of the results in [22| to a new class of 
stochastic differential equations. In |22, the author considers finite-dimensional realizations for HJM term 
structure models, see ^18) . This leads to the well known HJM drift condition for the underlying dynamics 
whereas we work with a different drift condition resulting into a different SPDE. 

The reminder of the paper is organized as follows. In Section 2 we set up a general framework for modeling 
of the futures curve. In Section 3 we analyze the case of finite-dimensional realization. The portfolio problem 
is recast as a finite-dimensional control problem in Section 4. 
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2. The Portfolio Problem 

2.1. Futures Price Dynamics. We start the analysis by defining tlie function space on wliicli tlie futures 
price curve t t-^ ft evolves. We use the Musiela parametrization and write the futures curve as a function of 
time-to-maturity y G R+. Following |Tl], see also j22i, we fix a parameter a > and denote as Ha the space 
of all absolutely continuous functions h : R+ — > R such that 

||/i|U:= (^\h{0)\^ + J\"y\h'iy)\^dyy < ^. 

Here, the derivative h' is understood in the weak sense. The space {Ha, \\ ■ \\a) is a separable Hilbert space 
for which the point evaluation h i— >■ 6y{h) : Ha — > R is a continuous linear functional for each y G R+ - see, 
e-g-, [22, Theorem 4.1. We observe from the definition of the norm || • \\a that the derivative of the functions 
in Ha decay at exponential rate whereas the actual futures price can very well be non-zero for large times to 
maturity. 

We denote as (S'f)f>o the semigroup of right shifts defined as Stf{y) = f{y + 1) for any given / G Ha- 
Furthermore, denote the differential operator ^ as A. Then we know that the semigroup (>5'f)f>o is strongly 
continuous with infinitesimal generator A, see [22], Thrm. 4.1. The domain of operator A is defined as 
V{A) ^{he Ha\ h' e Ha}- Furthermore, we define domains D(yl") {h G V{A''-^) \ A''-^h G V{A'''^)}, 
along with the intersection 'D{A°°) — Hi^i T^{A'^). Finally, an element h G 'D{A°°) is called quasi- exponential 
if the linear space spanned by the family {A"h}^^^ is finite dimensional. 

Let U he a. real separable Hilbert space. We assume that W be a Wiener process defined on a complete 
filtered probability space (51,7^, F,Q), where the filtration F = {J-^t}t>o satisfies the usual conditions, and 
taking values in U. Following the HJM-approach, we assume that Q is a (local) martingale measure. Denote 
the covariance operator of W as Q and the associated eigenvectors and -values as {ej}igi and {Xi}i^x, 
respectively. Here, I is a countable index set determined by the dimensionality of U. We assume that 
TrQ X] -^i < The family {ei}i^x is an orthonormal basis for the space U and we can express W as 

(2.1) Wt = ^VAlW^e„ 

where W are scalar Wiener processes. Furthermore, we denote Lf^ ~ Q^{U), where is the pseudo-square 
root of the covariance operator Q. To fix notation, we denote the space of all Hilbert-Schmidt operators from 
a given Hilbert space H^ to another Hilbert space iJ^ as Lhs(H^ , H^). 
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We assume that the risk- neutral futures price dynamics t ft are defined as the solution of the stochastic 
partial differential equation 



(2.2) dfti-) = AM-) dt + dWt = AM-) dt + Y. \^^<^^iM■)) dWl /o G 

where, for brevity, ai{h{-)) := E,(/i(-))(ei) G Ha- Here, the mapping S (Si)iei ■ -f^a ^ Lhs{U'^ ,Ha) 
satisfies appropriate measurability conditions, is Lipschitz and admits the uniform linear growth condition 
ll^(^)llLffs(c/'^ H ) — ^ ll^lla)' fo'" ^^ where C does not depend on h. These conditions guarantee the 
existence of a unique mild solution, see, e.g., [5], Chapter 7. The mild solution of (|2.2p can be expressed as 

M-)^Sth{-)+ f St-u {mui:))dWu) 
Jo 

= Stfoi-) + E ^ /* ('^^(/-(•))) ^W:- 

Furthermore, we assume that 



E 







< oo. 



This ensures that the mild solution is also a weak solution, see [15], Thrm 3.2. 

We need also a description of the futures dynamics under the market measure P. To this end we apply 
Girsanov's theorem, see, e.g., [9], Theorem 10.14. For a given ijj e U'^ , define the process VV as 

dWt = dWt - tpdt- 

Then W is a [/-valued Wiener process under an equivalent measure P defined via the Radon-Nikodym deriv- 
ative = £{—^ - W)t- Here, £ denotes the Doleans-Dade exponential, see [11], Chapter II. Furthermore, 
the covariance operator of W is Q. The futures price dynamics can be written under the measure P as 

dM-) = {AM-) + W) dt + dm 

(2.3) ■.= ,^{M-))dt + mt{-))dm- 

In financial terms, we can interpret the element ip as the market price of risk. Being an element of a Hilbert 
space, the market price of risk can have essentially richer structure than its finite-dimensional counterpart, 
a constant vector. Indeed, if the space U'^ is a space of functions of the time to maturity y, the market price 
of risk will also depend on time to maturity. Furthermore, we can interpret the term 'S{ft{-)){ip) as the risk 
premium. We observe immediately that the risk premium is proportional to the volatility of the futures price 
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curve. Our stochastic dynamics allow for a very flexible modelling of the risk premium, taking into account 
possible idiosyncraticites between futures contracts with different maturities. Since the P-dynamics of the 
futures price curve is completely determined by the volatility S and the market price of risk tp, we call the 
pair (S, V') a futures price model. 

2.2. Portfolios of Futures Contracts. As was mentioned already in the introduction, a futures contract 
is a derivative security written on the futures price ft- The initial market value of this contract is zero and 
it yields a cash flow during its lifetime generated by the fluctuations of the futures price. The proflts (losses) 
caused by the fluctuations are put in (drawn from) a margin account such that the brokers collateral remains 
on the required level over time. We assume that the investor can use the margin account also as a borrowing 
account with an interest rate same as the constant risk free rate of return. In order to describe the time 
evolution of the wealth generated by trading on futures contracts, consider flrst a single futures contract in 
discrete time. Denote the length of the time step as At. The parties enter a single futures contract at time t 
with the delivery at time T > t + At. This position is canceled at time t + At and a new contract is entered 
with the same time of delivery - this is the marking to market-pmcedme. The resulting profit/loss is the 
associated fluctuation of the futures price, i.e. 

ft+At{y - At) - ft{y) = ft+Atiy) - My) - Ut+At{y) - ft+At{y - At)), 

where y = T — t is the time to delivery at time t. When passing to the limit At — >■ 0, this gives the expression 

dftiy) ^ — dt. 

oy 

The wealth dynamics generated by a single futures contract investment can now be written as 
(2.4) dXt = (rXt - dt + dftiy), = 0. 



dy 

Here, r > denotes the risk free rate of return. It is worth pointing out that the wealth X can become 
negative. 

Consider now a general portfolio of futures contracts. A futures portfolio consists of future obligations 
to purchase or sell electricity for the futures price. As a starting point, we consider the portfolio strategies 
to be measure- valued processes 1 1->- r(t, ■). Here, the measure r(t, ■) gives the portfolio weights for contracts 
with times to maturity y > 0. In other words, it gives the net amount of obligations undertaken by the 
investor for any given time to maturity. Given a portfolio of futures contracts and the fact that it is costless 
to take positions in futures, the instantaneous payoff from holding this portfolio is determined solely by the 
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fluctuations of the associated futures prices. This resuhs into the weahh dynamics generated by an arbitrary 
portfoho strategy F which can be written as 



(2.5) dXf = (^rXf - ^ ^^r(t, dy)^ dt + dft{y)T{t, dy), = 0. 

To make the expression (12. 5p technically tangible, we assume that portfolios F : 1 1~> F(t, •) are F-progressively 
measurable processes taking values in the dual H* . Thus we can write the expression (|2.5p in the more compact 
form 

(2.6) dXf = (rXf - (m, ^^)) dt + (F(t),d/,(.)) , Xl - 0, 



dy 

where (•, •) is the bilinear pairing on H^. The wealth X^ can be written as 



(2.7) 



X,^ l^rXl -[T{s),-^Jjds + (F(.),C(-)) 

* (rXf + {T{s)Mftmm) ds + f (f(s), S(/,(.))dW.) , 







for all t e [0,r]. For the right hand side of (|2.7p to make sense, we additionally assume that the portfolio 
process F satisfies the L2(P)-condition 



(2., 



T 



F(i),S(/t(-))(^))dH +/ y^^,{V{t),a,{ft{-))fdt 



< oo. 



Portfolio processes F satisfying the condition (|2.8p will be called admissible. Given the time t and the value 
of the wealth process Xf = x, we denote the set of admissible portfolios as A{t, x). 

The portfolio strategies we are considering are not self-financing in the classical meaning of the word. 
Indeed, since it costs nothing to enter a futures position, the investor needs no initial wealth to set up a 
portfolio. However, the investor needs credit since the wealth (the amount of money in the external bank 
account, that is) can become negative. In this case, the investor needs to borrow money to balance the 
collateral, i.e. to infuse cash into the system. In this sense, the portfolio optimization problem considered 
here resembles more betting than traditional investment problem. The admissible portfolio strategies F 
describe the investor's allocation of wealth over the whole futures curve. This will include so-called roll-over 
strategies, where one invests in single a futures contract with a given time to maturity y. By picking more 
times to maturity yi, . . . , ?/„, one can create roll-over portfolios of futures contracts. 
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2.3. The Optimization Problem. The objective of the investor is to determine an admissible investment 
pohcy which maximizes the expected utihty at a terminal time T . To study this problem, we first propose 
an appropriate set of utility functions. We call a function m : R — ;> [— cxd, go) a utility Junction if 

(1) it is concave, nondecreasing and upper-semicontinuous, 

(2) the half line dom(u) := {a; £ R | u{x) > —00} is a nonempty subset of [0, 00), 

(3) u' is continuous, positive and strictly decreasing on the interior of dom(u) and UtHj-^qo u'{x) ~ 0, 

see [19) . p. 94. Given a utility function u, the portfolio optimization problem reads as 
(2.9) V{t,x)= sup i:,[u{X^))\Xf ^ x] , V{T,x) =u{x), 

with < i < T and a; > 0. 

3. FiNITE-DIMENSIONAL REALIZATIONS 

3.1. Invariant Foliation: A Characterization. In this section we assume that the real separable Hilbert 
space U is truncated to be finite-, say, n-dimensional. We can identify the truncated U with the euclidian 
space R". From a practical point of view, this corresponds, for example, to the case where the eigenvalues 
Xi are very small for all i > n. Thus the associated scalar processes W\ i > n contribute very little to the 
overall fiuctuations of the process W and, consequently, the processes Wi, i < n, can be identified as the 
principal components driving the futures curve ft, see, e.g., [7]. We denote the truncated (i.e., R"-valued) 
Wiener process as W. We remark that the market price of risk ip degenerates now into a constant vector 
(■01, . . . ,ipn) G R"- In |17) a two- factor model for the oil futures price dynamics is proposed, which corresponds 
in our context to assuming that the futures curve is driven by a two-dimensional Wiener process. On the 
other hand, in electricity there exist empirical and theoretical evidence for a high degree of idiosyncratic risk, 
which means that a high dimensional Brownian motion is required for the futures curve dynamics (see |20[ [T]). 

The first objective of this section is to pin down conditions on the volatility structure E under which the 
price dynamics given as the solution of the SPDE (|2.3p admit a finite-dimensional realization. The existence 
of a finite-dimensional realization means, roughly speaking, that for a given initial curve ho, there exists a 
nice family of manifolds in Ha such that the initial curve is on the initial manifold and that the evolution of 
the curve in confined to the family of manifolds. This problem setting has been studied extensively during 
the last decade or so, references include [31 IH [S] IHl HH HH [22] ■ In our analysis we follow the approach of [32]. 
We proceed by making the following definitions. 

Definition 3.1. Let V be a linear d- dimensional subspace of Ha- Then 
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(A) a family M := {A4t)t>o of linear submanifolds of Ha is called a foliation generated by V, if the there 
exist a continuously differentiable 6 : — > such that Ait = ^(^) + V, for all t > 0. Here, 9 is 
called the paranietrization o/M. Analogously, the tangent space is defined as TAit '■= 0' {t) + V for 
all t > 0, 

(B) a foliation M is called invariant for the SPDE (|2.3p if for all tg G R+ and Hq G Aito we have 
P(/t G Att„+t) = ^, t > 0, where ft is the (weak) solution for (|2.3p with /o — Hq. 

Definition 3.2. Let V be a linear d-dimensional subspace of Ha- Then the equation (j2.3p is said to have 
a finite dimensional realization generated by V if for all h G 'D(A) there exist a foliation M'' = {Alt)t>o 
generated by V with h G A^q which is invariant for (j2.3p . 

As we see from these definitions, the invariant foliation is the basic building block of a finite-dimensional 
realization. The next proposition gives necessary and sufficient conditions for invariance of a given foliation 
for the SPDE (|2.3p . This proposition is analogous to Theorem 5.3 in [22], where a similar characterization is 
proved when the underlying price dynamics are driven by a one-dimensional Wiener process. In our result, 
we consider the case where the driver is a multi-dimensional Wiener process. In addition, the drift term in 
the futures price dynamics (12. 3p cannot be handled immediately using existing results. Indeed, Theorem 5.3 
in |22| is concerned with the HJM approach of term structure modeling of forward rates which leads into the 
well-known HJM drift condition for the underlying dynamics, whereas we consider a different drift condition 
stemming from the fact that we analyze futures prices. For convenience, we recall the definition of z^: 

iy{h) = Ah + i:{h){iJ;), 

where h G VIA). Furthermore, we make the following assumption on the volatility structure. These as- 
sumptions allow us to use Theorem 2.11 in [22] (We point out that Theorem 2.11 in [22] is proved for a 
one-dimensional driving Brownian motion but we observe that the result holds also for a multi-dimensional 
Brownian driver by simply plugging it into the proof). 

Assumption 3.3. We assume that the components crj are continuously differentiable with aj{Ha) G {ft. G 
Ha' ■ h{oo) — 0}, for some a' > a, and that there exist L, A/ > such that 

\\a-j{hi) - crj{h2)\\a < L\\hi - h2\\a, hi, h2 G Ha, 

\\<Jjih)\\a<M, he Ha, 

for all j — I, ... , n. 
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Proposition 3.4. Let M be a foliation generated by the d-dimensional linear space V C Ha spanned by 
elements {vi}f^i. The following statements are equivalent: 



(A) The foliation M is invariant for the equation (|2.3 

(B) We have 



(3.1) 



(3.2) 



e{t) e V{A), for all t > 0, 
Vi e T>{A), for ffl/H = 1, . . . , d, 

and there exist functions (3 G C°^HR'+ x R.'^;R'') and k G C"^1(R+ x R'';R'* x R")stic/i that 

d 

V {B{t)) = e'{t) + J2 P^(t^ 0)t;, e TMu 

i=l 

(d \ d 

^(0 + ^^ZiVi = '^ViKij{t,z) e V, j = l,...,n, 
i=l ) 1=1 

for all {t, z) e R+ x R"^, and the elements Vk satisfy the ordinary differential equations 



(3.3) 



d d 



dy 



for all k — 1, . . . , c?. 



Proof. Necessity: Fix t > Q and assume that the foUation M is invariant for the equation (j2.3p . Then we 
know from (35], Thrm. 2.11, that M.t and, consequently, 6{t) and Vi are in T>{A) for all t > and i — 1, . . . ,d. 
Let h £ Mt and write h = 0{t) + J2i=i where z'* G R"*. A gain, we know from [22,, Thrm. 2.11, that 
i^{h) e TMt and (Tj{h) e V for aU j = 1, . . . , n. This yields 

d 

(3.4) a,(/.)=5]<^-('')^., 

where z^i^^^f g R'' for each j — 1, . . . ,ri. Define the linear isomorphism I : R'* ^> V as I{y) = ^i^j- 
Then we can write 

(3.5) K.j{t, z'^) := z'^^C'' = X-^ (^{t) + ^ zf 
for all J = 1, . . . , n. The claimed regularity properties of k follow from the assumptions on a. 
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To prove the claim on z/, we find similarly that 

d , d n d 



y 1=1 y i=l j=l i=l 



where z""^''^ G R'*. Then we find using (13.51) that 

(3.7) f3{t,z^) =Z-i ^±0{t) +J2 + E ^^'^^^ | -0'{t) 

Again, the claimed regularity properties of v follow from the assumptions on a. 

Finally, the differential equations p.3p follow by differentiating the expression p.6p with respect to 
for all fc = 1 , . . . , d. 

Sufficiency: Assume now that the conditions p.ip - p.3p hold and let h G A4t. Using [22], Thrm. 2.11, 
it suffices to show that 

Cd \ d 

e{t) + E ^^^^ - ^'(^) + E Z^'^^' ^'')^^ 
i=l I i=l 

to prove the claim. To this end, we observe first that (7j{0{t)) ~ ViKij{t, 0) for all j = 1, . . . , n. Thus we can 
write 

^(^w+E^N^) 

i=i y j=i 1=1 
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By using the differential equations p.3[) and changing the order of summation, we find that 



dy 



(i,0))V, 

d I " 



E /3fe(i,o)-E'^fe.(i.o)v, 



9Z; 

To proceed, we readily verify that 

n 

(3.9) =/3(i,2) - oi) I E«-j(i>^)V', 



On the other hand, since I is a linear isomorphism, the inverse X^^ : V — > R'' is also linear. Thus we observe 
from the expression p.9p that the partial derivatives z) — '^■ji'^T independent of z. 

Consequently, we find that 

f3{t, 0) - E 0)^^- + E 9^ ^) - E 

j=l i=l \ = 1 

is the Maclaurin series of the function /? — Y^^=i '^■j'^j- This proves the identity p.Sp . □ 



Proposition [33] gives a convenient characterization of invariance of a given foliation. First, it gives a set 
of ordinary differential equations p.3p which have the spanning functions Vi as their solutions. The coefhcients 
of these ODEs are characterized in terms of the original volatility structure E and the market price of risk 
ip. Furthermore, we observe that the volatility S must be of a very specific type in order to be associable to 
an affine realization. Indeed, each of the components aj must map the manifold A4 into the linear space V 
for all t. This means that the randomness in the price dynamics is confined to the finite-dimensional linear 
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structure V. Put differently, the price dynamics admitting a finite-dimensional realization can diffuse only in 
finitely many directions in the infinite-dimensional state space . 

We can identify the coordinate process Z driving the finite-dimensional realization using Proposition 
13.41 To this end, fix io G R-+ and h G M.to- Then we have a unique z G R'* such that h — 6{to) + X^iLi ^i'^i- 
Now, let Z be the strong solution of the Ito equation 

dZt = /3(to + t, Zt)dt + K{to + t, Zt)dWt, Zq = z, 

where /3 and k are given by Proposition [3]4l Then it is straightforward to verify using Ito's formula that the 
i/^-valued process / defined as 

d 

(3.10) ft := e(to + /o = 

is the strong solution of SPDE ((O)) . 



3.2. Interpretation of the Coordinate Process. The coordinate process Z has a priori no intrinsic 
economical interpretation. However, we can equip it with one. To this end, we observe that the futures curve 
can be decomposed as 

ft = TTvft + TTyi ft = TTv/t -|- 'rT^±Mt+to ■ 

Here, the latter equality (which holds almost surely) follows from the invariance of the foliation M and the 
fact that the projection TTy±A4t+to is a singleton set. We point out that the projection ny±A4t+to can be 
used as a parametrization of the foliation and that it is the unique parametrization which is in for all 
t > 0. According to p.lOp . we can rewrite the futures price as 

d 

(3.11) ft^ny^Mt+to+Y^Zlv,. 

1=1 

Denote a basis of the subspace V"*^ as {wj}°^i and let A : Ha R'' be a linear continuous operator such 
that A(V) = R''. Then we can write 

Cd \ oc d 

^ Zlv. = ^ c>A{w,) + J2 ZlAM 
i=l / j=l i=l 

(3.12) 

OO d d d / oo \ 

= E ^' E ^^'^^(^») + E ^tHv^) = E + E ^'^^^ ^(^^)- 
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Define the matrix Q as Qij — Ai{vj), wliere i,j = 1, . . . , d. Then we can write 

oo 

z, = A{f,)e-'-Y,cib,.. 

In other words, we observe that the state process Z can be associated to the image of / in arbitrary linear 
operator A modulo an affine transformation depending of A and the invariant foliation M. This means, in 
particular, that the coordinate process Z becomes an observable quantity as a functional transformation of 
quantities which are observed from the futures price curve. 

Example 3.5. Consider the operator A given by benchmark contracts for times to maturity i/i, i = 1, . . . , d, 
defined as 

Here, 5 is the evaluation functional on Ha, = h{y). Then we can write 

d oo 

j=i i=i 

for all i — l,...,d. The dimension d gives the number of benchmark contracts needed to reconstruct the 
whole futures curve t ft via the identity p.lOp . Indeed, we need d benchmark contracts in order to have a 
connection between the coordinate process Z and the benchmark contracts f-iUi)- This connection is given by 
the invertible d x d-matrix Q with coefficients Ai(vj) ~ Vj{yi), i,j — \,...,d. 

3.3. Existence of Finite-dimensional Realization: Sufficient Conditions. In the previous subsection, 
we proved a characterization of the invariance of a given foliation M for the SPDE (|2.3p . In this subsection 
we use this result to give sufficient conditions for the existence of a finite-dimensional realization. To this 
end, consider the futures price dynamics given by (j2.3p when the volatility is of the form 

p 

(3.13) aj{h) = he Ha, 

i=l 

for a\\ j = 1, . . . ,n. Here, the functions Vi, i = 1, . . . ,p, are linearly independent and each $ij maps Ha into 
R. We assume that the functionals $ij are twice continuously differentiable and that there exist constants 
L,M > such that 

||%(/ll) - $y(/l2)IU < -^11^1 - ^2|U, hi,h2 G Ha, 

<M, he Ha, 
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for alH = 1, . . . ,d and j = 1, . . . ,n. This guarantees that we can use Proposition l3.4l The equation (|2.3p can 
now be written as 

(n p \ n p 

j=l i=l J i=l i=l 

where are scalar P-Brownian motions. We remark that the expression p.lSp gives actually a necessary 
condition for the existence of a finite-dimensional realization, see |22) . Lemma 3.2. Thus this form of volatility 
is the most general we can consider in this framework. The following proposition, which is analogous to 
Proposition 6.2 in [22], gives sufficient conditions for the existence of a finite-dimensional realization. 

Proposition 3.6. Assume that the volatility structure of (|2.3I) is given by (|3.13p and that the functions Vi, 
i = 1, . . . ,p are quasi-exponential. Then the equation (j2.3p admits a finite- dimensional realization. 

Proof. Since the elements Vi, i — 1, . . . ,p are quasi-exponential, the linear space 



r:=0Span|— z;, : n>ojcV{A) 



is, by definition, finite-dimensional. Furthermore, we observe that 

(3.15) -^v&Y, veY. 

ay 

Set d := dim Y and choose Wp+i , . . . ,Vd Y such that {vi , . . . , w^} is a basis of Y . 
Fix an arbitrary Hq G 'D{A). First, define the function 9{t) : R+ Ha as 

(3.16) 6{t) := Stho G P(A). 
and the function k e C°-\R+ x R''; R'^ x R") as 



(3.17) Kij{t,z) 



0, i — p -\- 1, . . . ,d. 



for all J = 1, . . . , n. We observe from (|3.17p that the latter condition in (|3.2p is satisfied. Furthermore, define 
the function ^ e C°'\Tl+ x R'';R^) as 

71 d 

(3.18) I3t{t,z) = ^Kij(t, z)'0j -I- y^gjfcZfc, 
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where chosen, due to (I3.15p . such that 

d 

(3.19) —v, = ^vjaji, 

for alH = 1, . . . , d. With this specification, we verify by differentiating expression (I3.18P with respect to Zk 
that the elements vt satisfy the differential equations (|3.3I) . Finally, the definition of (3 coupled with the fact 
that -^0{t) = 9'{t) = ha [y + 1) implies that the former condition in (I3.ip is also satisfied. □ 



Proposition l3.6l and its proof give us not only sufficient conditions for the existence of a finite-dimensional 
realization but also a recipe for its construction. To illustrate this, we consider the example from Section 6 
in [3]. 

Example 3.7. Assume the driving Brownian motion in (j3.14p is one- dimensional, the volatility functional 
$ = 1 and the quasi- exponential function Vi{x) = xe^°'^ for a > 0, i.e., that a = ui = vi. In addition, let 
for simplicity — 0. We observe that the space V = Span{ui,U2}, where V2{y) — e^"^, so the dimension of 
the affine realization is 2. It is a matter of differentiation to show (see p.l9|) ) 



-a 
1 



Thus the state variable Z is given, due to (I3.17P and (j3.18p . as the strong solution of the two-dimensional ltd 
equation 

dZt = AZtdt-hBdWt, 

where B = (1,0)^ and W is a C^-Brownian motion. Now, given the initial ho G T>{A), the affine realization 
of futures price reads as 

ft{-) = hoi- +to+t) + Zlvi{-) + Zfv2{-). 
It is worth pointing out that the process Z becomes an Ornstein- Uhlenbeck process with solution given by 

Zt = exp{At)Zo + [ exp(A(t - s))B dW^ . 
Jo 

As the eigenvalue of the matrix A is —a, that is, negative, Zt becomes stationary. Many popular .spot and 
futures price dynamics of commodities and energy have stationarity as a crucial property, see [2 for examples 
in energy. The model in the current example gives rise to a hump-shaped futures price curve, which is relevant 
in oil markets, say (see |16] ). 



16 
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4. Portfolio Optimization Revisited: Finite-dimensional State Variable 

In this section we recast the optimization problem (j2.9l) into a finite-dimensional setting using the 
identity (13.101) . Assume that the initial futures price curve can be expressed as /o — St^ho + X^iLi -^j'^^i where 
ho G 'D{A) and to G R+. Then the affine futures price curve dynamics is given by 

d 



ft = StStgho + ^ Zlvi. 



Recall that the wealth X is given by the expression (|2.7p . To rewrite this for the affine realization, we first 
observe that 

dh 

dft = StSto -7T-dt + ^ dZlv^. 



Then we can rewrite the wealth as 



(4.1) Xf = 1^ [rX^ 
It is a matter of differentiation to show that 



c, c- dfs _ dho dho ^ dv, _ A dv. 

By coupling this with the expression (I3.19p . we observe that the equation (|4.ip can be written as 

Xf = rxUs + (r(.), g dziv, - 1: Zl^ds^ 

= [ rXlds+ [ /r(s),^dZ>,-^Z:^i;,a,,dA 
•^0 -^0 \ ^=l j=l I 

rXlds + {t{s),Y1^v, \dZl-Y,a,,Zlds 

rt d / d \ 

rXUs + J^I1 7t \ dZl - J2 a^JZids \ , 

where the real valued processes 7* := {T,Vi) satisfy the condition (|2.8p . Finally, by using identity (13. 18^ . we 
find 

„t d n 



(4.2) 



(4.3) 



•^0 -^0 ^=l j = l ^ 

= f rX]ds+ f -fs{dZs-AZ,ds), 
Jo Jo 
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where is a P- Wiener process and the matrix A is given by (|3.19l) . By plugging this into the utility 
maximization problem (|2.9p . we reduce the initial infinite-dimensional control problem into a classical finite- 
dimensional control problem. We write this problem in the Markovian form 



(4.4) 



V{t, X, z) = supE [m(XJ)|X7 = x,Zt = z\. 



We can use classical control theoretic techniques to study this problem, see, e.g., [14]. To illustrate this, we 
derive the HJB-equation for this problem. Assuming that the function V is sufficiently smooth, Ito's formula 
yields 



(4.5) 

dV{t,X'<,Zt)^ 



d 71 



«>i=i 



k=l 



dt 



Vt{t, A7, Zt) + V,{t, Xl Zt) [rx + Y."^,Y. + II ^t, Zt)l3,{t, Zt) 

\ i=l fe=l / i=l 

^ d n 1 ^ ^ 

+ -Vxx(t,X^,Zt)'^'^JiKij{t,Zt)Kji{t,Zt)ji + - ^ VzjZ, {t, A7, Zt) Kjkjt, Zt)Kki {t, Zt) 

i=l k=l 

d n 
+ ^ Vxz^{t,X^, Zt) ^ Ktk{t, Zt)Kki{t, Zt)ji 
i=l k=l 

d n d 

+ vAt, A7, Zt) Y^^Y. ^^k(t, Zt)dWt + Y ^t, Zt) Y ^^k{t, Zt)dWt. 

i=l k=l i=l k=l 

By a standard martingale argument from stochastic control, see, e.g., [2], this yields the HJB-equation 

{/ d n \ d 

Vt{t,x,z) + Vx{t,x, z) rx -f ^7i^ z)?/;*; J + ^ V^, (i, x, z)/3i(t, z) 

V 1=1 k=l ) i=l 

^ d n 1 ^ 

(4.6) + -VUt,x,z)YT. 

{t,z)Kji{t,z)ji+ - ^ Vz^zj {t, X, Z) Kjkjt, z)Kkr{t, z) 



i=l k=l 



k=l 



'^V^zA'tyXy z)'Y l^ik{t, z)Kki{t, z)ji I = 0, 



k=l 



where V{T, x, z) = u{x) for all a;, z G R. As a consequence, we can formulate the following verification result 



Proposition 4.1. Assume that we have 
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(1) a smooth function V : [0, T] x R*^ x R — > R with V{T, z, x) = u{x) such that 

Cd n \ d 

i=l fe=l / i=l 

^ d n 1 ^ ^ 

+ ■^Va:x{t,X, z)^^'JiKij{t, z)Kji{t, z)ji + 7^^ Vz^z,{t, X, z)^ K-jkit, z)Kki(t, z) 



i=l k=l 



k=l 



'^VxZi{t,X, Z)'^ Kikit, z)Kki{t, Z)ji < 0, 



i=l /c=l 

/or t, 7, X, z, 
(2) a function 7* : [0, T] x R x R'' ^ R'^ swc/i that 

Cd n \ d 

rx + ^ 7* , X, z) ^ Kik {t,z)^k\ + V;. {t,x,z)fii{t,z) 
i=l k=l 
d n 



4=1 



i't,zh*{t,x,z) 



i=l k=l 



i,i=l 



fc=l 



+ '^Vxz,it,X, z)'^Kik{t, z)Kkt{t, z)j*{t,X, Z) = 0, 
i=l /c=l 

for all t, z, X. 

(3) a wealth process X* corresponding to an admissible control 7 such that 7f = 7*(i, , z) for all t > 
and z e R''. 

Then the function V solves the control problem (14. 4p (i.e. V ~ V) and an optimal control is given by 7. 

Proof. Let 7 be an admissible control and X'* the associated wealth process. Then Ito's formula yields 
(4.7) 

dV{t,XlZt)^ 



k=l 



V i=l / i=l 

d n 1 ^ ^ 

+ Y ^t''' H '^^^^^^ Zt)Kk^{t, Zth^{t) 

i=l k=l 

d n d n 



dt 



1=1 fc=i 



4=1 



fc=l 
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Since the term on the right hand side inside the angled brackets is negative, we obtain 



= V{T, XJ, Zt) < V{t, X], Zt 



V i=i k=i 

d n \ 

+ Y,V^As,X2,Z,)Y,'^^k{s,Z,)dW^ . 



i=l k=l J 

Since 7 is admissible, the stochastic integral is a true martingale and, consequently, E[m(XJ)|X7 — x,Zt = 
z] < V{t,Xf , Zt). This shows that V dominates the solution V. To show that V < V and that 7 gives an 
optimal control, we use the same argument coupled with the conditions (2) and (3). □ 



Proposition 14. II gives a set of sufficient conditions for a given smooth function to coincide with the value 
of the control problem (|4.4p . Moreover, it gives an identification of an optimal control. Even though the 
control 7 does not appear to have an direct economical interpretation, it can be traced back to, for example, 
to benchmark futures prices following Example 13.51 In this example, the d-dimensional coordinate process 
Z was identified with an affine transform of d benchmark futures prices. On the other hand, as we observe 
from the equation (|4.3I) . the process 7 can be understood as "portfolio weights" on d-dimensional diffusion 
dynamics t i~> k(s, Zs) (^dWs + ipds^ . These dynamics can extracted from the evolution of the coordinate 
process Z. 

We end our paper with discussing a futures price dynamics typical for energy markets. Suppose that 
the dynamics is given as 

(4.8) dft{x) = {Aftix) + aiV-i + a2e-""V2) dt + ai dWl + aae-'^" dW^ , 

for a two-dimensional Brownian motion W and cri,(T2,a positive constants. For simplicity, we also assume 
the market price of risk i = 1, 2 to be constants. In this model, we identify $12 = $21 = 0, $11 = cri and 
$22 = cr2- Furthermore, vi{x) = 1 and V2{x) — exp(— aa;), which obviously are quasi-exponential functions 
spanning a linear space of dimension 2. As v'i{x) = and v'2{x) — -~av2{x), we find the dynamics of the 
two-dimensional process Z to be 

dZt = + AZt )dt + BdWt, 

with = ((TiV'i, o'2'/'2)' and 










0-1 





A = 




B = 






-a 







0-2 
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We conclude by the analysis above on finite dimensional realizations that 
(4.9) ft{x) = ho{x + to + t) + Z\ + e-'^'-'Zl, 

for some to > 0. By splitting into the two components of Z, we find that is a drifted Brownian motion 

dZ^ =0-1 V-i dt + (Ji dWl, 

and Zf is a mean-reverting process, 

dZf = (cr2V'2 - aZf) dt + (72 dWf. 

Optimizing a portfolio invested in the futures price curve will with this model be equivalent to optimizing a 
portfolio investment in two "assets" with dynamics Z. As exp(— ax) tends to zero when x tends to infinity, 
we find that ff {x) ho{x + Iq + 1) + Z} for large values of x. Hence, an investment in Z^ can be viewed as 
holding a portfolio position in a futures with long time to maturity, that is, a position in a contract in the 
far end of the futures curve. The "asset" Z^ can then be interpreted as the difference between a futures far 
out on the curve (being Z^) and one with short time to maturity. Our investment problem will therefore be 
to select optimally a portfolio of contracts in the short and long end of the curve. 

The model in (j4.8|) . or (j4.9p . can be viewed as the implied futures price dynamics from a two- factor spot 
model. In fact, following [T7], we can assume that the spot price of some commodity is given by 

St = m + z\ + zi, 

where A(t) is some deterministic seasonality function. This spot price model will be an arithmetic analogue 
of the dynamics proposed for oil by [17 . Here, Z^ is interpreted as the short term variations of the oil spot 
price, while Z^, the non-stationary part, is the long-term trends in oil prices including inflation and extinction 
of reserves. This corresponds to the view of investing in the long and short end of the futures curve. 
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